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(57) ABSTRACT 

A system that simultaneously measures the translational tem- 
perature, bulk velocity, and density in gases by collecting, 
referencing, and analyzing nanosecond time-scale Rayleigh 
scattered light from molecules is described. A narrow-band 
pulsed laser source is used to probe two largely separated 
measurement locations, one of which is used for reference. 
The elastically scattered photons containing information 
from both measurement locations are collected at the same 
time and analyzed spectrally using a planar Fabry-Perot inter- 
ferometer. A practical means of referencing the measurement 
of velocity using the laser frequency, and the density and 
temperature using the information from the reference mea- 
surement location maintained at constant properties is pro- 
vided. 
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SPATIALLY-AND TEMPORALLY-RESOLVED 
MULTI-PARAMETER INTERFEROMETRIC 
RAYLEIGH SCATTERING SYSTEM AND 
METHOD 

5 

RELATED APPLICATIONS 

The present application is a continuation of PCT Applica- 
tion No. PCT/US1 1/035580, filed May 6, 201 1, which claims 
priority to U.S. Provisional Application Nos. 61/424,502, 10 
filed Dec. 17, 2010, 61/346,599, filed May 20, 2010, and 
61/332,511, filed May 7, 2010, the contents of which are 
incorporated herein by reference. 

STATEMENT REGARDING FEDERALLY 15 

SPONSORED RESEARCH OR DEVELOPMENT 

This invention was made with government support under 
contracts NASW-99027, NNX07AC3 2 A/2007 and 

NNX08AB31 A/2008 awarded by NASA. The government 20 
has certain rights in the invention. 

FIELD OF THE INVENTION 

The present invention relates to spatially and temporally 25 
resolved multi-parameter Interferometric Rayleigh scatter- 
ing. 

BACKGROUND OF THE INVENTION 

30 

Molecular-based diagnostics techniques capable of obtain- 
ing simultaneous multiple fluid properties such as tempera- 
ture, density and velocity are critically important for charac- 
terizing the flows within airbreathing engines, such as 
scramjet engines and scramjet-turbine combined cycle 35 
engines for hypervelocity vehicles. At the other end of the 
velocity range, such techniques are needed to help in the 
development of the next generation of subsonic, low noise 
aircraft. Statistical correlations between these properties will 
lead to a more detailed understanding of the complex flow 40 
behavior as well as aid in the development of new turbulence 
models for scramjet engine flow path designs, and more accu- 
rate aircraft noise prediction tools. Required are instanta- 
neous and simultaneous measurements of temperature, 
velocity, density, pressure, and chemical composition at mul- 45 
tiple points in a gas when the spatial (hundreds of microns or 
less) and temporal (hundreds of nanoseconds or less) scales 
of the turbulent fluctuations are resolved. Specifically of 
interest are the mean properties and their turbulent variances 
and co-variances as well as the probability density llinctions. 50 
To date, diagnostics techniques to meet all these requirements 
on a single platform do not exist, but many linear and non- 
linear optical techniques are available for partial achievement 
of this goal. One of these linear techniques, the Interferomet- 
ric Rayleigh Scattering technique (IRS) (discussed in Biv- 55 
olaru, D., Danehy, R M., Gaffney, Jr. R. L., and Cutler, A. D., 
“Direct- View Multi-Point Two-Component Interferometric 
Rayleigh Scattering Velocimeter,” AIAA-2008-0236, 46th 
Aerospace Sciences Meeting, Reno, Nev., Jan. 9-12, 2008), 
has been used for low-as well as for high-speed supersonic/ 60 
hypersonic non-reacting and combusting flows yielding tem- 
porally and spatially resolved simultaneous measurements of 
two-components of bulk velocity. 

Recent efforts have attempted to obtain the ro-vibrational 
temperature, and the species content of N2 and 02 from 65 
measurements with CARS (coherent anti-Stokes Raman 
spectroscopy), simultaneously with velocity with IRS. Those 
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efforts include Danehy, P. M., Magnotti, G., Bivolaru, D., 
Tedder, S., and Cutler, A. D., “Simultaneous Temperature and 
Velocity Measurements in a Large-scale, Supersonic, Heated 
Jet,” Paper 1193, 55th JANNAF Propulsion Meeting, Boston, 
Mass., May 12-16, 2008. Also, Bivolaru, D., Lee, J. W., Jones, 
S. B., Tedder, S., Danehy, P. M., Weikl, M. C., Magnotti, G., 
and Cutler, A. D., “Mobile Rayleigh — CARS Instrument for 
Simultaneous Spectroscopic Measurement of Multiple Prop- 
erties in Gaseous Flows,” 22nd International Congress on 
Instrumentation in Aerospace Simulations Facilities 
(ICIASF), Monterey, Calif., June, 2007; and Bivolaru, D, 
Danehy, P. M., Grinstead, K. D., Jr., Tedder, S., and Cutler, A. 
D., “Simultaneous CARS and Interferometric Rayleigh Scat- 
tering,” AIAA-2006-2968, 25th AIAA Aerodynamic Mea- 
surement Technology and Ground Testing Conference, San 
Francisco, Calif., Jun. 5-8, 2006. However, this work needs 
improvements in instrument precision, as well as additional 
simultaneous measurements, the gas density and a measure- 
ment of translational temperature from IRS. 

SUMMARY OF THE INVENTION 

Accordingly, it is an object of the invention to provide a 
multi-parameter IRS technique for simultaneously obtaining 
the translational temperature, the bulk velocity, and the den- 
sity in gases at elevated temperatures in a non-intrusive man- 
ner. In one non-limiting illustrative embodiment, the inven- 
tion measures these properties in a H,-air combustion flame 
(Hencken burner), relative to simultaneous measurements in 
a reference N 2 -filled cell. The measurements are compared 
with the adiabatic flame theory predictions and the applica- 
bility of this method to turbulent flows is discussed. 

In accordance with these and other objectives, a system is 
provided which simultaneously measures the translational 
temperature, bulk velocity, and density in gases by collecting, 
referencing, and analyzing nanosecond time-scale Rayleigh 
scattered light from molecules. A narrow-band pulsed laser 
source is used to probe two largely separated measurement 
locations, one of which is used for reference. The elastically 
scattered photons containing information from both measure- 
ment locations are collected at the same time and analyzed 
spectrally using a planar Fabry-Perot interferometer. A prac- 
tical means of referencing the measurement of velocity using 
the laser frequency, and the density and temperature using the 
information from the reference measurement location main- 
tained at constant properties is described. A simplified Gaus- 
sian distribution model to the scattered light spectra is used to 
obtain the flame properties. Corrections to this model are 
applied at lower gas temperatures when the simplified Gaus- 
sian approximation is no longer suitable. The near-zero mea- 
sured velocity as a function of the measured flame tempera- 
ture, and a comparison of the measured flame density and 
temperature with the perfect gas law are presented. 

These and other objects of the invention, as well as many of 
the intended advantages thereof, will become more readily 
apparent when reference is made to the following description, 
taken in conjunction with the accompanying drawings. 

BRIEF DESCRIPTION OF THE FIGURES 

FIG. I shows the instrument configuration for simulta- 
neous multi-point measurement of velocity, the gas transla- 
tional temperature, and density at two largely-separated spa- 
tial locations. Legend: L-lens; M-planar mirror; BD-beam 
dump; PBS-polarization dependent beam splitter/combiner; 
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DP-Dove prism; F-Pi-Fabry-Perot interferometer; em-CCD- 
electron avalanche multiplication CCD camera; mv-mea- 
surement volume. 

FIG. 2 shows Fabry-Perot interferograms containing mul- 
tiple signals of elastically scattered light from two largely 
separated measurement locations, one of which is used for 
reference. The bottom horizontal pattern (expl) is the signal 
obtained at the location of a Fl 2 -air Flencken burner flame: (a), 
in air at 293 K outside the flame, and (b), in the flame at about 
1560 K. The top horizontal pattern (exp2) is the signal from 
the secondary measurement location situated in a low-pres- 
sure N, gas cell at room temperature. This signal is used to 
reference the measurement of temperature and density per- 
formed at the flame location. The concentric ring pattern 
(refl) is the spectra of unshifted laser light used to reference 
the Doppler frequency of all signals and facilitate computa- 
tion of the gas bulk velocity at both locations. 

FIG. 3 shows the combined Rayleigh and laser spectra 
obtained, (a), in atmospheric air at -293 K, and (b), in FI 2 -air 
flame at 4>~0.5 (1637 K). The intensity scale is shown in 
arbitrary units (a.u.). Data contain ten instantaneous measure- 
ments to emphasize the trend (solid circles). Only one dataset 
(empty circles) is fit with a combined Gaussian-Lorenzian 
function (solid black line). The residual of the fit is shown 
with a line W. 

FIG. 4 is a chart showing normalized spectral properties 
obtained simultaneously in H2-air Hencken burner flame 
(<|)~1.0, and <|>~0.5), and in the reference N2-cell relative to 
atmospheric properties: (a), the spectral bandwidth, (b), the 
area of the spectra, and (c), the Doppler shift frequency. The 
inset on the top of the figure is a top view of the central section 
of the burner. 

FIG. 5(a) is the measured gas density in air and in flame as 
a function of the measured temperature, compared with the 
perfect gas law (solid curve). FIG. 5(b) is the measured bulk 
velocity in air and in flame as a function of the measured 
temperature. The gas density and temperature are normalized 
with respect to the ambient air density (1 .2 kg/m3) and tem- 
perature (293.5 K), respectively. The velocity is normalized 
by the maximum measurable velocity (3 km/s). Data obtained 
at atmospheric pressure during a full scan through a H 2 -air 
Hencken burner flame at about 4>=0.5 . The burner stabiliza- 
tion co-flow is nitrogen gas. 

FIG. 6 is a combined single-shot interferogram of signal 
(horizontal pattern), frequency reference (concentric ring 
pattern), and density reference from an optical fiber (semi- 
circular pattern). 

FIG. 7 is an instantaneous image containing two interfero- 
grams of different spectral resolutions. 

DETAILED DESCRIPTION OF THE PREFERRED 
EMBODIMENTS 

In describing the preferred embodiments of the present 
invention illustrated in the drawings, specific terminology is 
resorted to for the sake of clarity. However, the present inven- 
tion is not intended to be limited to the specific terms so 
selected, and it is to be understood that each specific tenn 
includes all technical equivalents that operate in a similar 
manner to accomplish a similar purpose. 

The interferometric Rayleigh scattering (IRS) technique 
uses a narrow-band polarized laser source of wavelength k 0 
for probing gas molecules situated at a given measurement 
volume through photon scattering. Both elastic Rayleigh and 
inelastic Raman scattering of photons occurs. The inelastic 
scattering arises from changes in the rotational and vibra- 
tional states of the molecules during the scattering process. 
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This process is weak and occurs in a different wavelength 
band than the laser wavelength. The Rayleigh scattering pro- 
cess is directly related to the spatio-temporal motion of the 
gas molecules and does not change the internal eneigy of the 
5 molecules. It is orders of magnitude stronger and occurs near 
the laser wavelength. The number of photons scattered are 
proportional to the incident laser energy. Receiving optics 
directly collect the scattered light (near kg) in a small solid 
angle, optically compare it with laser light for frequency 
to reference, and analyze it spectrally with a high-resolution 
Fabry-Perot interferometer and a CCD camera detector. 

The details of the spectral structure are related to the state 
of collisions between gas molecules traditionally divided in 
two regimes: hydrodynamic or collision dominated, and 
15 kinetic or collisionless regime (Rayleigh regime), without a 
clear border between them. The shape of the spectra is related 
to a non-dimensional ordering parameter y=a/Ka, where 
a=p/p is the collision frequency (p is the gas pressure and p is 
the shear viscosity) and Ka is the acoustic spatial frequency 
20 (K is the magnitude of the resultant wave vector K and a is the 
speed of sound). For “s” type scattering, the angle between 
the electric field vector E of the incident light of wave vector 
k 0 (vertically polarized) and the direction of the scattered 
light of wave vector k^ (horizontal direction) is p=90°. The 
25 scattering angle 0 is the angle between the direction of k 0 
vector and the direction of k 5 vector along the optical axis of 
the receiving optics. 

With these notations the resultant wave vector is K=k 5 -k 0 . 
Assuming that the wave vectors magnitudes are about the 
30 same Ik^Nk,-,!, the magnitude of the resultant wave vector is 

K=2k 0 sin(0/2) (1) 

where k 0 =lk 0 l=2 n/k 0 . For very low density gases y is negli- 
gibly small and the Rayleigh spectra is approximated with a 
35 Gaussian shaped spectra. For moderate density gases, when 
collisions between molecules start to dominate (y— 1), the 
spectra include a central peak (Rayleigh) and two Lorentzian- 
shaped sidebands (acoustic scattering or Brillouin). The Ray- 
leigh contribution arrives from thermal fluctuations dissipat- 
40 ing by thermal conduction. For large ordering parameters, the 
line width in frequency is given by D r K 2 , where D r is the 
thermal diffiisivity. The Brillouin sidebands are symmetri- 
cally Doppler shifted about the central Rayleigh peak due to 
acoustic scattering from small pressure fluctuations in the gas 
45 (propagating at the speed of sound). They are located at 
±Ka/2ji, and their line width is given by TK 2 , where T is the 
sound attenuation coefficient. Then, the ratio of Rayleigh to 
Brillouin integrated intensities is (y-1), where y is the ratio of 
specific heats. This approximation fails in the kinetic regime 
50 when y is negligibly small. 

At elevated temperatures for atmospheric pressure (or low 
pressures at room temperature) (y«l), if we neglect the 
acoustic Brillouin scattering, the velocity distribution of mol- 
ecules is Maxwellian. For one velocity component this is 
55 given by 

yKM^r^expt-K-v) 2 /,^ 2 ) (2) 

where ^ is the velocity of a molecule, v is the bulk velocity of 
molecules, 

60 

Wxr (2k B T/m) m (3) 

is the full width at 1/e height of distribution f(^), T is the gas 
temperature, k s is the Boltzmann constant, and m is the 
molecular mass. The bulk velocity v is what is observed at 
65 macroscopic scale. For a gas at rest the bulk velocity is zero 
and the average molecular random velocity u=(^-v) due to 
deviations of ^ from v is also zero. 
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Photons of frequency u from the incident probing beam are 
scattered from molecules moving in all directions with the 
random velocity u, and are shifted in frequency according to 
the Doppler effect 

dv (K9 u)/2ji (4) 5 

Therefore, for a component of velocity along the resultant 
wave vector K one can relate the velocity distribution of the 
molecules to the scattered light frequency by rewriting Eq. (2) 

as to 

f(v)=(7iw} 1/2 exp (-(dv/w) 2 ), (5) 

where the new spectral width in frequency is 

w=(2k B T/m) 112 sin(0/2 )/X 0 . (6) 

In a gas mixture (as for example in air or in a flame) the 15 
spectral distribution of the scattered light is a sum of distri- 
bution with contributions from each molecular species i 
weighted by their molar fraction and the differential scat- 
tering cross-section a,. Therefore one can write 


Hi (7) 

l(v) = CIoN^XiO-iMv), 


where I 0 is the laser intensity, and N is the total gas number 
density. The constant of proportionality C depends on the 
laser polarization, and wavelength; the size of the probe vol- 
ume; and the solid angle of the light collection optics (col- 
lecting power), the collection direction, and the detector 30 
quantum efficiency. 

To determine the gas properties, the scattered light spectra 
and the laser spectra are analyzed with theoretical models by 
least-square fitting techniques using the laser spectra as ref- 
erence for Doppler shift calculations. Research efforts to 35 
extract the properties of gases by analyzing the spectra of 
scattered light with Fabry-Perot interferometers include in 
diluted gases, in low-speed combustion flames, and in high- 
speed rocket external flows. The use of kinetic models to 
obtain spectral information requires prior knowledge of gas 40 
parameters (that are not measured by the instrument) such as 
the internal specific heat per molecule, the ratio of shear to 
bulk viscosity, and the ratio of shear viscosity to thermal 
conductivity. This model is difficult to use since these param- 
eters are not precisely known and there is limited data avail- 45 
able for the species and range of temperatures required for 
measurements. In this work a simplified Gaussian-Lorentzian 
analytical model is used, easy to apply to the data, instead of 
a more detailed computational Rayleigh-Brillouin kinetic 
model. Then, corrections of these measured quantities via 50 
calibrations are applied to obtain quantitative information 
within an acceptable margin of error as it will be shown later. 

FIG. 1 is a non-limiting embodiment of the instalment 
design in accordance with a preferred embodiment of the 
invention. The system 100 references all measurable proper- 55 
ties, the velocity, the temperature, and the gas density. This is 
performed with an optical setup that simultaneously records 
the spectral properties of Rayleigh scattered light from two 
largely separated measurement locations, one of which is 
used for reference. To increase the incident laser energy prob- 60 
ing the measurement volume, and implicitly to increase the 
scattered signal, this system uses for probing two collinear 
laser beams slightly delayed in time for probing the gas 
sample. The laser source is a seeded dual Nd: YAG laser, with 
a wavelength, X 0 , of 532 nm, variable timing between pulses, 65 
10 Hz repetition rate, nine nanoseconds pulse width and 
0.005 cm -1 bandwidth. 
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As shown, the probing laser beam (of wave vector k 0 ) 102 
is focused at the measurement volume (mvl) 106 by the h 1 
lens 1 04 and collected and re-collimated for further use by the 
L 2 lens 108. The L7 lens refocuses the recovered laser beam 
after passing through mvl to generate the second measure- 
ment volume mv2. Thereafter, the remaining beam energy is 
lost in the beam dump BD. The elastically scattered light 
pulse (expl) (of wave vector kj 120 is collected and colli- 
mated from the measurement volume mvl 106 by the L 3 lens 
112 (which can have an f/4, 250-mm focal length, and is 
located 250 nun from the measurement volume mvl 106). 
The beam diameter of this signal is reduced to match the 
analyzing optics aperture diameter (-15 mm) using the beam 
reducer/expander formed by lenses L 4 114 and L 5 116, and 
directed using the mirror Mj 118 toward the polarization 
dependent beam combiner PBS 122. The beam combiner 
PBS 122 facilitates the optical mixing of the signal with the 
un-scattered light (refl) 124 used for frequency referencing. 

In the present configuration, the polarization dependent 
beam combiner PBS 122 is used to also combine this beam 
(expl ) 120 with a second beam (exp2) 130 collected by the L s 
lens 136 from a secondary measurement volume mv2 140. 
This measurement volume 140 is located in a gas cell with 
known gas properties that is used forproperties referencing as 
will be shown later. Similarly, tins beam (exp2) 130 is reduced 
and collimated by the lenses L g 134 and L 10 132, and directed 
by mirror M 3 126 toward the combiner PBS 122. The light 
polarization of both signals is adjusted (not shown) such that 
the PBS 122 reflects the second beam exp2 130 together with 
the transmitted main beam 120 toward the Fabry-Perot inter- 
ferometer (F-Pi) 1 42. The reference laser light 124 needed for 
frequency referencing is optically mixed with the signals at 
the PBS 122 by passing un-shifted laser light through the 
mirror M 3 126 (or any other mirror, such as mirror M x 118. 
Therefore the interferogram formed by the F-Pi interferom- 
eter 142 and imaged by the L 6 lens 144 on the CCD camera 
146 contains the spectra of expl 120, exp2 130, and refl 124. 

FIG. 2 shows examples of interferograms containing mul- 
tiple Rayleigh signals (expl and exp2) obtained in the cell 
(top horizontal pattern) and at the flame location (bottom 
horizontal pattern), in room air (293 K in FIG. 2(a)) before 
flame ignition, and in an H 2 -air flame at equivalence ratio 
(j>=0.5(1560 K in FIG. 2(b)), respectively. The circular pattern 
is the interferogram of the laser light used for frequency 
reference (refl). To avoid the eventual overlap of the closely 
spaced spectra expl and exp2 the pattern of exp2 could also 
be rotated in the image plane (for example by 90° or any other 
angle) using the Dove prism DP 128. 

As mentioned earlier, in this arrangement the energy 
deposited at the measurement volumes is increased by spa- 
tially overlapping the two collinear laser beams slightly 
delayed in time. Additionally optical aberrations are intro- 
duced in the optical path of the focusing lens through lens 
tilting. This beam shaping is achieved by tilting the focusing 
lens Lj slightly, 5° to 22° depending on the laser energy, and 
thereby introducing an astigmatism aberration of the beams at 
their foci. The tilting is performed about an axis perpendicu- 
lar to the plane containing the incident laser beam and the 
signal collection direction. The elliptical beam allows more 
energy to be deposited in the flow without gas breakdown. 
This results in larger Rayleigh signal intensity if collected 
from a direction in the plane of the beam and the maj or axis of 
the ellipse, at the cost of reduced spatial resolution in the other 
dimension. In the experiments reported here a total of 385 mJ 
in two laser pulses separated in time by about 200 nanosec- 
onds and a ratio of major to minor axis of the elliptical beam 
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of about 10 at its first focus were used. The measurement 
volume is about 0.05 mmx0.2 mmx0.5 mm. 

Both L 3 and L 8 lenses 112, 136 collect scattered light from 
the measurement locations at an angle 0=90° with respect to 
the corresponding probing laser beam direction. The Fabry- 5 
Perot interferometer 142 is made of solid glass with a free 
spectral range (FSR) of 8 GHz and a total finesse of about 22. 
This defines the maximum measurable Doppler shift and 
therefore the dynamic range for measuring velocities to about 
±1.5 km/s when no prior velocity direction is known, or to 10 
about 3 km/s if the velocity direction is known. The range of 
measurable temperatures and densities in gases is from cryo- 
genic temperatures to about 2500 K (limited at high tempera- 
tures by the overlap of consecutive orders of interference), 15 
and from 0.1 kg/m 3 to about 2 kg/m 3 (limited by the CCD 
dynamic range), respectively. 

A non-limiting example of the invention is discussed to 
illustrate the invention. The first measurement location is 
situated in a non-premixed near-adiabatic H 2 -air Hencken 20 
burner flame at atmospheric pressure (also used to calibrate 
the CARS measurements). The burner is placed vertically at 
the first measurement volume (mvl) in a direction perpen- 
dicular to the scattering plane. The hydrogen fuel and the air 
oxidizer are assumed to mix immediately above the burner 25 
surface. The burner surface is a 25.4 mmx25.4 nun rectangu- 
lar section of hypodermic needle tubes for fuel and honey- 
comb passages for oxidant surrounded by a rectangular 
region of co-flow honeycomb design (a total of 36.5 
mmx36.5 nun square surface). The PI, and air equivalence 30 
ratio 4> is varied from 0. 1 to 1 , by varying the flow rates (up to 
0.5 1/sec). A constant flow of N, of 0.4 1/sec was directed 
through the co-flow nozzle at all fuel-air mixtures to improve 
flame stability and provide an almost dust-free environment 
for the instrument. The molecular content in the co-flow- 35 
flame interface is considered unknown. The flame tempera- 
ture and species molar fractions in the flame region were 
calculated from the adiabatic flame theory. 

The primary species are N,, O,, H,0, and H,, with the 
minor species (Ar, CO, OH, etc.) being neglected. The main 40 
product of combustion is water (H,0) is in the gaseous state 
with negligible contribution from condensate. The perfect- 
gas law for the mixture is p/p=(R/m)T, where m is the 
molecular weight of the mixture (kg/kmole), p is the gas 
pressure (N/m 2 ), p is the density (kg/m 3 ), and T is the gas 45 
temperature (K), and the perfect-gas constant R is 8314.5 
J/(kmole)(K). The experiments are performed in room air 
(m al ,.=28.97 kg/kmole) at the atmospheric pressure p ajr 
(101324.72 N/m 2 ), temperature T al> (293.5 K), and density 
p atr (1.2 kg/m 3 ). The axial velocity of the flow (vertical direc- 50 
tion) was estimated from flow rates at room temperature to 
about 10 m/sec, but higher velocities of this order of magni- 
tude are expected during combustion. The radial velocity 
component (horizontal direction, measured in this experi- 
ment) is considered to be much smaller than the axial velocity 55 
(near zero). 

The second measurement volume is situated in a constant 
properties glass cell containing N, at room temperature 
(293.5 K). The glass cell is connected with manual valves and 
plastic tubing to a feedstock of N, gas and to a vacuum pump . 60 
The initial air in the cell is removed with the vacuum pump 
and the N, is let in to the desired pressure before measure- 
ments. A pres sure ofp 0 =26.5 kPa (198.5 torr) was used for the 
reported experiments. Lowerpressure was used to ensure that 
the density and the y-parameter are similar to those in the 65 
combustion flow, which made the signals from both channels 
comparable. 
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Verification of the instrument is performed by obtaining all 
scattered light spectral properties simultaneously in the atmo- 
spheric flame and nitrogen (N,) cell . Single-shot spectra 1 (u) 
obtained in the air and H,-air flame are shown in FIG. 3. The 
data are shown with solid and empty circles. The best fit to the 
data shown with empty circles using the sum of a Gaussian 
function for the scattered spectra and a Lorentzian function 
for the reference laser spectra is shown with a continuous line. 
Ten single-shot spectra are overlapped to the fitted data (solid 
circles) to emphasize the trend. 

FIG. 3 (a), and FIG. 3(b), show the spectra obtained in 
atmospheric air at the conditions described earlier, and in the 
flame at (j>=0.5 (-1637 K), respectively. The energy of the 
laser was maintained approximately constant during these 
measurements though any variations in laser intensity would 
be corrected by the data obtained in the reference cell. The 
parameter y for air at atmospheric conditions is 0.8 so the 
spectra in FIG. 3(a) should theoretically contain Brillouin 
sidebands Doppler shifted by about ±0.9 GHz around the 
Rayleigh center peak. The residual between the experimental 
data and the Gaussian fit is shown in the plot at the bottom of 
the figure (curve designed W). This indicates that the Gaus- 
sian fit function over-predicts the wings and under-predicts 
the top half of the distribution. Therefore the distinctive 
maxima of this residual function approximately detect the 
Brillouin peaks locations and the integral of its absolute value 
approximately shows the acoustic scattering contribution to 
the spectra. A theoretical analysis of this integral function of 
the y-parameter, shows that this contribution to the measured 
Rayleigh spectra is percentwise similar with the one found 
experimentally by Cattolica et al. (“The Interpretation of 
Spectral Structure of Rayleigh Scattered Light from Combus- 
tion Gases,” AIAA-1 976-0031, 14'* Aerospace Sciences 
Meeting, Washington, D.C., January 267-28, 1976) for dif- 
ferent gases. This relation is used to correct the measured 
spectra for the Brillouin contribution at lower gas tempera- 
tures. 

For a Gaussian profile, the gas translational temperature is 
calculated from the spectral width. Normalizing Eq. (3) with 
respect to the known values of properties, the non-dimen- 
sional gas translational temperature becomes T m =W m 2 rn, 
where W m is the measured non-dimensional width of the 
spectral profile, and in the molecular weight of the gas, which 
is a function of the gas composition, normalized by the 
molecular weight of air. This relation holds at moderate to 
high temperatures, but generates large errors without know- 
ing the gas composition, and when the spectral profile devi- 
ates from a Gaussian function as described earlier. Since the 
molecular weight is not measured, it was estimated by com- 
paring the measured width with the computed spectral width, 
which is a function of the molecular weight of the gas (and 
temperature). Then, the computed molecular weight of the 
gas is used as an approximation of the molecular weight of the 
gas existing at the time of measurement. In doing so it is 
assumed that for certain measured spectral widths only a 
small range of molecular weights are possible. For example, 
in the calculation of T m , from the spectra shown in FIG. 3, the 
broadest Rayleigh spectra shown in FIG. 3(b) cannot be 
attributed to cold gases (higher molecular weights) that will 
produce the spectra shown in FIG. 3(a). The maximum error 
possible in these computations without considering this 
approximation of the molecular weight is about 16% in the 
whole range of measurable temperatures presented here. The 
actual error of the measurement obtained by estimating 
molecular weights as explained in this approximation is con- 
sidered to be only a fraction of this maximum error. 
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To obtain the bulk velocity of the gas, the spectra are 
analyzed for the average spectral shift frequency relative to 
the laser spectra. From Eq. 4, if all molecules move with the 
bulk velocity v, the spectra is Doppler shifted from the laser 
frequency u 0 by (u- , u o )/D o =2(v/c)sin(0/2), where u 0 (=c/X 0 ) 5 
is the frequency reference of the laser. 

In a usual configuration for measuring gas density, to keep 
the measurement unaffected by the fluctuation in laser energy 
from shot to shot, a reference measurement of the laser energy 
is implemented by recording a part of the signal or the laser to 
light. The solution for a reference adopted in this work is a 
secondary measurement volume located in a glass cell con- 
taining a gas with known pressure, temperature, and compo- 
sition. By normalization with the cell signal total intensity A 0 
the laser energy fluctuations from pulse to pulse during the 15 
measurement are minimized. Additionally, the measured cell 
parameters provide simultaneously a reference for known gas 
density, temperature, and velocity (zero for stagnant gas). 
Integrating over the solid angle of the collecting optics and 
rewriting in dimensionless quantities, Eq. (7) becomes 20 
p„,=C 2 (A m /A 0 )/(cr/m m ), where A m A 0 is the non-dimensional 
integral of the signal, p m is the gas density normalized by p air , 
and o=E J N '(x i o j )/a ajr is the non-dimensional scattering cross 
section of the gas weighted by the molar fraction of species x,- 
The constant C 2 is found from simultaneous measurements of 25 
known gas properties in the cell (N 2 ) and in air at the mea- 
surement volume. 

Tire total measured signal is computed by integrating the 
spectral fit function (solid line in FIG. 3 minus the reference 
Lorentzian function). For y«l, the integral can be approxi- 30 
mated by (assuming a Gaussian function) the maximum spec- 
tral intensity multiplied by full-width half maximum and 

(2jt) 1/2 . 

FIG. 4 shows the profiles of measured quantities in the 
Hencken burner flame at equivalence ratio of <j>= 1 .0 (line X) 35 
and 4>=0.5 (line Y), relative to the known quantities in the cell 
(Z trace). Data is extracted from single-shot interferograms 
using Gaussian and Lorentzian fit functions as approximation 
for the Rayleigh and laser spectra. FIG. 4(a) shows the nor- 
malized spectral bandwidth W m measured at the flame loca- 40 
tion. The normalization factor is the spectral width W 0 
obtained at this location in air at room temperature (-293.5 
K). The width obtained in the N 2 cell, about 9% narrower than 
the width measured in the air, can also be used for reference. 
FIG. 4(A) shows the integral of the signal intensity measured 45 
at the flame location, relative to the integral of the signal 
obtained in the cell, corrected for the cell signal intensity 
relative to the flame signal intensity. The core flame cross 
section, where the flame properties are almost uniform, is 
about 20 mmx20 mm. The measurement volume was scanned 50 
along a line from one side to the other over the center of the 
burner (shown in the inset image of FIG. 4) in a plane situated 
at about 38.5 mm above the burner surface (similar measure- 
ments of rotational-vibrational temperature using CARS 
were perfonned by Hancock (“Nitrogen and Hydrogen 55 
CARS Temperature Measurements in a Hydrogen/ Air Flame 
Using a Near- Adiabatic Flat-Flame Burner,” Combustion and 
Flame, Vol. 109, No. 3, 1997, pp. 323-331)). The Doppler 
shift frequency obtained in flame at equivalence ratio of 4>=0.5 
relative to the maximum measurable frequency (defined by 60 
the interferometer) is shown in FIG. 4(c). The simultaneous 
properties measured in the reference cell are plotted with red 
symbols. The normalization factors for the cell properties in 
FIG. 4(a) and FIG. 4(A) are the measured cell properties, the 
spectral width and the integral of the signal in the cell at one 65 
atmosphere, respectively. The spatially resolved measure- 
ments of the relative spectral width show a slight flame asym- 
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metry in this property as the measurement volume is scanned 
across the flame region (FIG. 4(a)). This is evidenced more at 
higher flame temperatures (blue) than at intermediate ones 
although in both cases they correlate approximately with the 
position of the burner nozzles revealing the discrete (flamelet) 
structure of the flame. In the core region of the flame the 
measured average temperature is 2323 K±41 5 K (±one stan- 
dard deviation) on the left side and 2083 K 343 K on the right 
side of the axis of flame at (|>=1 .0. These average values are 
about 2.2% and 12.4%, respectively, lower than predicted. 
The large standard deviation of the measurement is attributed 
partly to the variation with position (during the scan) in flame 
temperature and composition due to the flamelet structure and 
partly to random measurement error. At moderate tempera- 
tures ( — <))=0.5), when the flame is more uniform, this variation 
is smaller than 80 K. Because the measurements were not 
perfonned at a fixed location in the flame, the instrument 
random error in measuring temperature is not directly quan- 
tifiable at this time. 

The measured density in air, and flame at t()=l are 1.23 
kg/m 3 +0.05 kg/m 3 and 0.13 kg/m 3 ±0.007 kg/m 3 , respec- 
tively. These mean values are higher than predicted, by about 
2.5% in room air and by about 4.5% at the flame. The higher- 
than-expected value of density in the flame suggests also a 
flame of lower temperature than predicted as well as mea- 
sured from the spectral width. 

The gas density, and bulk velocity computed from spectral 
properties obtained in air and the flame (data shown in FIG. 4) 
during the full scan of the flame at <)>=0.5, are shown in FIG. 
5 as a function of the computed temperature. The gas density 
and temperature are normalized with respect to the ambient 
air density (1.2 kg/m 3 ) and temperature (293.5 K), respec- 
tively. The velocity is normalized by the maximum measur- 
able velocity (3 km/s) defined by the optical setup, i.e. the 
interferometer free spectral range. 

As shown in FIG. 5(a), the measured gas density depen- 
dence on the measured translational temperature (blue sym- 
bols), agrees very well with the perfect gas law (solid curve) 
computed for this case. In this computation the gas tempera- 
ture and composition is assumed to be only a function of <|) so 
the gas molecular weight can be inferred as a function of the 
temperature outside of the flame core where <|) is unknown. 
The measured velocity dependence on the flame temperature 
is shown in FIG. 5(A). The rms (root mean square) variations 
that include true innate instrument errors and spatial and 
temporal fluctuations of the flow. As expected, the mean 
velocity is near zero, and the mis fluctuations (±one standard 
deviation) of velocity increases with increasing temperature 
from about 1 1 m/s at the room temperature (and in the cell) to 
about 27 m/s in the flame core (and increases to 39 m/s at 
(j>= 1 ). It should be noted that this instrument was not intended 
to measure near- zero velocities, but rather supersonic/hyper- 
sonic ones (up to 3 km/s). Therefore the error of measuring 
near zero velocity is less than 1 % of the dynamic range of the 
instrument. To accurately measure low velocity flows, a 
higher resolution interferometer (lower free spectral range) 
and a backward scattering direction (more sensitive to Dop- 
pler shift) could be selected. 

To our knowledge, these high resolution temporal and 
spectral measurements of simultaneous gas translational tem- 
perature, velocity, and density using interferometric Rayleigh 
scattering are the first to be performed in H 2 -air combustion 
flames. 

Thus, the invention provides a non-intrusive optical diag- 
nostic technique for temporally and spatially resolved mea- 
surements of multiple properties in reacting and non-reacting 
gases. The system employs a narrow -band pulsed laser source 
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at 532 nni for simultaneous probing of molecules through 
elastic light scattering at two largely separated measurement 
locations, one of which is used for reference. The collected 
photons containing information from both measurement 
locations are analyzed spectrally with a planar Fabry-Perot 5 
interferometer. The reference measurement is located in a 
glass cell containing N 2 gas at low pressure (26.5 kPa) and 
room temperature (293.5 K). The measured cell spectral 
properties provide simultaneously a reference for known gas 
density (the spectral intensity), translational temperature (the to 
spectral width), and supplementary near-zero gas velocity 
(zero average Doppler shift). Additionally, the cell spectral 
intensity provides the normalization factor for the laser 
energy fluctuations critical for the measurement of density. 
These spectral properties are used to calculate, by nonnaliza- 15 
tion, the gas density and translational temperature at the mea- 
surement location. The range of measurable temperatures and 
densities in gases is from cryogenic temperatures to about 
2500 K (limited by the interferometer), and from 0. 1 kg/m 3 to 
about 2 kg/m 3 (limited by the CCD dynamic range), respec- 20 
tively. The maximum dynamic range for measuring velocities 
is about ±1.5 km/s when no prior velocity direction is known, 
or about 3 km/s if the velocity direction is predictable. 

In a laboratory demonstration, the measurement location is 
situated in a non-premixed near-adiabatic H 2 -air Hencken 25 
burner flame at atmospheric pressure. The referenced mea- 
surements at different flame conditions compare favorably 
with the flame temperature and the gas density (up to the 
stoichiometric flame temperature of 2377 K) computed with 
the adiabatic flame theory. The random errors of measuring 30 
the flame properties, the translational temperature, velocity, 
and the gas density at the maximum measurable gas tempera- 
ture, are estimated to less than about 17%, 39 m/s and 5%, 
respectively. These errors include the spatial and temporal 
variation of these properties across the flame. 35 

The invention demonstrates simultaneous referencing 
(calibration) while performing measurements of velocity, 
density, and temperature with the interferometric Rayleigh 
scattering technique described here. In particular, the inven- 
tion can be used for simultaneous multi-property turbulence 40 
studies of subsonic, supersonic, and hypersonic, nonreacting 
flows and predictable composition reacting flows in the range 
of temperatures, turbulence levels, and errors described here. 

In another embodiment of the invention, an optical setup is 
implemented in which only one Rayleigh scattered light pulse 45 
(signal expl) is collected from the measurement volume 
(mvl) by the LI lens, and collimated and directed toward the 
PBS. Referring to FIG. 1, the PBS is used to inject the beam 
into a recycling optical path formed by the input mirror of the 
first interferometer F-Pi situated on the first optical axis and 50 
the input mirror of a second interferometer (E2 — not shown) 
situated on a second optical axis perpendicular to the first 
optical axis. This interferometer replaces mirror M 3 and is 
adjusted such to reflect the signal received from PBS. A 
quarter-wave plate (QWP — not shown) located between the 55 
PBS and F-Pi rotates the polarization of the input beam expl 
to circular polarization. 

The signal pulse polarization from mvl is adjusted such 
that the PBS sends it toward the QWP and interferometer F-Pi 
first, and reflects it towards the QWP, PBS and the second 60 
interferometer (E2) afterwards. A small fraction of this light 
is transmitted through either interferometer, but most of it 
(depending on interferometer reflectivity) is reflected back 
into the opposite interferometer. In the usual arrangement, not 
containing the second interferometer this rejected light from 65 
F-Pi is lost. In the disclosed optical arrangement the reflected 
light by the F-Pi input mirror is reflected back on the optical 
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axis again by the interferometer input mirror of the second 
interferometer. At the second reflection the QWP rotates the 
polarization to the initial input polarization and the remaining 
beam energy is lost on the first optical axis. This is demon- 
strated in Bivolara, et al., “Multipoint Interferometric Ray- 
leigh Scattering using Light Recirculation,” AIAA-2008- 
3708, 26th AIAAAMT-GT Conference, Seattle, Wash., Jun. 
23-26, 2008, using a mirror instead of a second interferometer 
(see also Bivolaru et al., U.S. Pat. No. 7,414,708). Optical 
means (not shown) directs the light passing through E2 
toward L6 to generate a second interferogram. 

In a further aspect, the interferometer F-Pi is selected for 
the maximum range of measurable velocities while the sec- 
ond one for a much lower range (or vice versa), e.g., with the 
free spectral range of 8 GHz (low spectral resolution) and 
3.75 GHz (high spectral resolution), respectively. The inter- 
ferometer F-PI is used to remove the ambiguity in the sign of 
the velocity being measured, while the second one measure it 
accurately. At higher temperatures this interferometer cannot 
be used due to the overlap of the consecutive interference 
orders, the second interferometer suffices. Similarly, at very 
low gas temperatures (<290 K) when the Rayleigh spectra are 
narrow, the second interferometer suffices. The dynamic 
range and accuracy of the instrument can be adjusted by 
selecting the free spectral ranges of the interferometers, 
depending on the application. 

In a further aspect, for the disclosed implementation a 
second direct -viewing collecting optic (similar to lenses LI, 
L2 and L3) collects and collimates Rayleigh light (ref2) from 
a second measuring volume situated in a controlled environ- 
ment, FIG. 1 . The light polarization is adjusted and directed 
such that the PBS reflects it toward F-Pi. After one pass 
through the interferometer the reflected signal cannot be 
recovered and is lost. To avoid the overlap of spectra, the 
reference pattern of ref2 is rotated 90° in the image plane 
using the Dove prism DP. Therefore the interferograms gen- 
erated by both interferometers, contain the spectra of expl, 
refl (probing laser spectra), and ref2 (reference spectra), and 
the spectra of exp2 (the high resolution Rayleigh spectra) and 
refl, respectively. 

In another embodiment of the invention, a single direct- 
viewing optic is used to collect light from two (or more) 
largely-separated measurement volumes. One (or more) mir- 
rors) or equivalent optics placed in close proximity to the 
lens LI location, receives and transmits (toward the lens) light 
from additional measurement volumes. One measurement 
volume is used for direct probing as in the new or prior art 
while the secondary measurement volumes are used for mea- 
surements of density, temperature, and velocity in conjunc- 
tion with one or more reference cells, or known gas flows. An 
example of a density (velocity and temperature) reference cell 
layout is shown in FIG. 1; or by imaging optical fiber(s) that 
collect signal(s) from widely separated measurement or ref- 
erence location(s). An example of an interferogram used to 
combine two signals, one directly from a measurement vol- 
ume (vertical pattern), and one from an imaged optical fiber 
used for density reference (three semi-circular sectors) is 
shown in FIG. 6 . The optical fiber collects light remotely from 
a density reference cell. The circular pattern is the interfero- 
gram of the laser light used for frequency reference. 

FIG. 6 is for a combined single-shot interferogram of sig- 
nal. frequency reference, and density reference from an opti- 
cal fiber. The single shot image shows an interferogram of a 
signal collected from the measurement volume, an interfero- 
gram of signal from a secondary measurement volume used 
as density reference, and an interferogram of the probing laser 
light used as a frequency reference. Thus, a single-shot inter- 
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ferogram is provided which combines a signal from a direct 
view optic as described previously, a frequency reference, and 
a density reference from an optical fiber. 

In another embodiment of the invention shown in FIG. 7, 
the output fringe pattern of both interferometers are optically 5 
combined in a single beam and imaged onto a single high- 
sensitivity CCD camera. FIG. 7 shows the instantaneous 
image containing two interferograms of different spectral 
resolutions. A first interferogram 220 is shown having a free 
spectral range (FSR) of 3.75 GFlz (F-Pi). A second interfero- 10 
gram 224 is shown having an FSR of 8 GFlz (E2). The hori- 
zontal pattern 222 shows the additional non-spectrally 
resolved density signal delivered parallel with expl from the 
measurement volume. The vertical pattern 228 shows the I5 
spectrally resolved signal from the cell for referencing the 
density measurement (ref2). The frequency reference pat- 
terns (refl) 226 are also shown. This provides better accuracy 
at lower velocities while preserving a large dynamic range. 
And, eliminates the Doppler shift sign ambiguity at higher 20 
velocities. 

In another embodiment of the invention, the signal is fur- 
ther increased by increasing the laser energy deposited at the 
measurement volume without inducing the gas breakdown. A 
laser sheet or an elliptical beam at its focus allow more energy 25 
to be deposited in the flow, resulting in larger Rayleigh signal 
intensity if viewed from a parallel direction and more accu- 
rate measurements, at the cost of reduced spatial resolution in 
the other dimension. This can be achieved by tilting the focus- 
ing lens L l slightly (5-30 degrees) thereby introducing astig- ,(1 
matism, by using a combination of cylindrical and spherical 
lenses, or by a pair of anamorphic prisms or similar to gen- 
erate a laser sheet. 

In another embodiment of the instrument a laser sheet , . 
generated by tilted lens Lj is viewed from a perpendicular 
direction such to permit two-dimensional measurements of 
parameters without decreasing the equivalent signal obtained 
from a single circular cross-section beam. 
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The description and drawings of the present invention pro- 
vided in the paper should be considered as illustrative only of 
the principles of the invention. The invention may be config- 
ured in a variety of size, shapes, and optical elements, and is 
not intended to be limited by the preferred embodiment. 15 
Numerous applications of the invention will readily occur to 
those skilled in the art. Therefore, it is not desired to limit the 
invention to the specific examples disclosed or the exact 
construction and operation shown and described. Rather, all 
suitable modifications and equivalents may be resorted to, 20 
falling within the scope of the invention. 

The invention claimed is: 

1. A system for analyzing properties of gaseous media, 
comprising: 

a laser source providing a narrow-band pulsed laser which 25 
probes a location in the gaseous media and a reference 
location; 

a first optical element collecting scattered light from the 
gaseous media; 

a second optical element collecting scattered reference 30 
light from the reference location; and, 

a beam combiner optically mixing the scattered light col- 
lected by said first optical element with an unscattered 
light reference and the scattered reference light by said 
second optical element. 35 

2. The system of claim 1, further comprising an interfer- 
ometer receiving the mixed laser from said beam combiner. 

3. The system of claim 2, further comprising a CCD camera 
for imaging the light from said interferometer. 

4. The system of claim 3, further comprising a processor 40 
for analyzing the light from said interferometer to identify 
properties of the gaseous media. 

5. The system of claim 4, wherein the system operates in 
real time to provide a diagnostic of gaseous media flows. 

6. The system of claim 1, wherein the system operates in a 4 - 
range of temperatures, velocities, and gaseous media densi- 
ties. 

7. The system of claim 1, wherein the gaseous media is in 
the form of clusters of molecules (Mie scattering). 

8. A system for spatially and temporally resolving multiple 50 
parameters of a gaseous media sample, the system compris- 
ing: 
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a probing laser source providing a narrow-band pulsed 
laser source; 

a first optical element for focusing the laser on the gaseous 
sample; 

a second optical element collecting and collimating scat- 
tered light from a gaseous sample; 
a first beam reducer/expander reducing a diameter of the 
collected and collimated scattered light from the gas- 
eous sample; 

a polarization dependent beam combiner; 
a first mirror directing the scattered light from said beam 
reducer/expander to said beam combiner; 
a third optical element collecting and re-collimating the 
laser passing through the gaseous sample; 
a second mirror directing the laser from said third optical 
element; 

a fourth optical element focusing the recovered and colli- 
mated laser on a reference measurement volume with 
known gaseous media properties; 
a second beam reducer/expander reducing the diameter of 
the scattered light received from the second optical ele- 
ment; 

a third mirror directing the scattered reference light from 
said second beam reducer/expander to said beam com- 
biner, whereby said beam combiner optically mixes the 
scattered light from said first mirror with an unscattered 
light used for frequency referencing, and with the scat- 
tered reference light from said third mirror; and, 
a first Fabry -Perot interferometer receiving the mixed light 
signals from said beam combiner. 

9. The system of claim 8, wherein the unscattered reference 
light passes through said third mirror toward said beam com- 
biner. 

10. The system of claim 8, wherein said first, second, third, 
and fourth optical elements each comprise a lens or a combi- 
nation of lenses. 

11. The system of claim 8, wherein said first and second 
beam reducer/expanders each comprise two lenses. 

12. The system of claim 8, wherein the said system probes 
one location in a gaseous media using simultaneously two 
interferometers in a signal recycling configuration, the signal 
recycling configuration comprising: 

wherein said third optical element comprises a second 
Fabry-Perot interferometer receiving the light signals 
from the said beam combiner, wherein the second Fabry- 
Perot interferometer is adjusted perpendicular to an opti- 
cal axis; and, 

a quarter-wave plate located between the beam combiner 
and the first Fabry-Perot interferometer, wherein the 
said quarter-wave plate modifies linear polarization of 
an incoming light beam to circular polarization. 





